Can the frequency dependent isobaric specific heat be 
measured by thermal effusion methods? 
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Abstract. It has recently been shown that plane-plate heat effusion methods devised for wide-frequency specific-heat 
spectroscopy do not give the isobaric specific heat, but rather the so-called longitudinal specific heat. Here it is shown that 
heat effusion in a spherical symmetric geometry also involves the longitudinal specific heat. 
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INTRODUCTION 

The frequency-dependent specific heat is one of the most 
fundamental thermoviscoelastic response functions char- 
acterizing relaxation of liquid structure in highly viscous 
liquids. The measurement of this quantity as an alterna- 
tive to enthalpy relaxation studies was conceived more 
than two decades ago |[ll|2l- 

When the specific heat is frequency-dependent one 
faces experimentally the problem of separating out the 
trivial frequency dependence from the slow propagation 
of heat. This can be solved in two ways. One can go 
to the thermally thin limit [2^ ^ where the sample is 
small compared to the thermal diffusion length lo- Id is 
inversely proportional to the square root of the frequency 
and typically 0.1mm even at IHz. So this condition can 
be difficult to fulfill over a wide frequency range unless 
the sample is very small jit]. 

The other - effusion - approach is to choose the sam- 
ple size much larger than the thermal wavelength. This 
is the thermally thick limit and it is easier to realize over 
a wide frequency range |4]. Both methods have to take 
due account of the stresses coming from the material sup- 
porting the sample associated with different thermal ex- 
pansion coefficients. The effusion method suffers addi- 
tionally from thermal stresses produced within the liquid 
itself. It is only the latter problem we consider in this pa- 
per. 

In the effusion methods one typically produces a har- 
monically varying heat current Re{fft)e"*"} at a surface 
in contact with the liquid. The corresponding tempera- 
ture response RelroDg""'} on the very same surface is 
measured. Since the response is linear in the stimulus it is 
convenient to introduce the complex thermal impedance 
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The thermal impedance of a sample of volume V and 
volume specific heat c is 
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in the thermally thin limit. 

If on the other hand planar thermal waves effuses from 
a plate of area A into a liquid the thermal impedance is 
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in the thermally thick limit fll]. Here X is the thermal 
conductivity. 

It has allways tacitly been assumed that measurements 
done at ambient pressure are isobaric and that the c of 
formula ([3j is Cp. However it was formerly stated JH] and 
recently shown that the ordinary heat diffusion with a 
complex diffusion constant does not describe the experi- 
mental situation adequately. For unidirectional heat effu- 
sion it was shown that the effective specific heat mea- 
sured is the so-called longitudinal specific heat c/(a)) 
which is between the isochoric, cv(ft)) and isobaric, 
Cp{(X)) specific heats. Denoting the adiabatic and isother- 
mal bulk moduli by Ks{co) and Kt{co) respectively and 
the shear modulus by G{co) one can write the adiabatic 
and isothermal longitudinal moduli asMs{(o) = Ks{(o) + 
4/3G(«) and Mrico) =Kt{(o)+4/3G{(o). Now Clio) 
is related to cv{co) as 0] 
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whereas Cp{co) is related to cv{co) as 
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In an easily flowing liquid G{co) is negligible since l/co 
is large compared to the Maxwell relaxation time, Tm, 
and there is no difference between c/ and Cp. However in 
a highly viscous liquid near the dynamic glass transition 
the shear modulus becomes comparable to the bulk mod- 
uli and the diffence between c/((b) and Cp{co) becomes 
significant. 

We show below that the same is true for heat effusion 
in a spherical symmetric geometry. Here one also obtains 
ci{co). Furthermore in spherical geometry one can also 
get the heat conductivity A and thus get c/(a)) absolutely. 
The planar unidirectional method in fact gives only the 
effusivity, y/Xci{co); that is ci{co) is determined only to 
within a proportionality constant. 



THERMAL AND MECHANICAL 
COUPLING 

The general equations 

One cannot treat the diffusion of heat independently 
of the associated creation of strains or stresses. Let the 
temperature field, T{r,t) be described in terms of the 
small deviation dT{r,t) — T{r,t) — To from a reference 
temperature Tq and denote the displacement field by 
u = u(r,r). Dealing with relaxation is most conveniently 
done in the frequency domain. Thus time dependence of 
the fields is given by the factor e*', s = ico. Considering 
only cases where inertia can be neglected the equations 
that couple temperature and displacement are 

MrV(V • u) - /3v/V5r - GV X (Vxu)=0 (6) 
cvsdT + l5vTosV-u-XV^5T = 0. (7) 

Here the isochoric pressure coefficient )3v(a)) is defined 
in the constitutive equation for the trace of the stress 
tensor a 

ltr{a)=KTVu-Pv5T. (8) 



The isobaric case 

If the trace of the stress tensor is constant in time 
then the term jSyTo^V • u in equation (|7| becomes 
Tq^y/Kt5T. Since TqPy/Kt = Cp — cy equation (|7]i now 
becomes the ordinary heat diffusion equation 
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decoupled from the displacement field and with a diffu- 
sion constant involving the isobaric specific heat 



Cn 



It is usually assumed that thermal experiments on liq- 
uids with a completely or partially free surface will be at 
isobaric conditions. However this is only true if the shear 
modulus G can be neglected compared to buUanodulus 
Kj. This condition fails near the glass transition and the 
full coupled problem of equations (|6]l and (jTjl has to be 
considered. 



Radial heat effusion from a spherical 
surface into an infinite media 

We would like to show here that the inherent problem 
of measuring Cp is not only confined to one-dimensional 
heat flow in the geometry considered in 1^ where the 
associated displacement field is forced to be longitudinal. 

The longitudinal specific heat also emerges in the ther- 
mal impedance against effusion out from a sphere. In the 
spherically symmetric case V x u vanishes. If we denote 
differentiation with respect to r by a prime (|6]l and (|7| 
becomes 
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Mrir^ir^u)')' ^l5v5T' = (11) 
8T + Tolivsr-^{r^uy - Xr-^{r^5T')' ^0. (12) 



Define now the longitudinal specific heat. 



the heat diffusion constant. 



and the wave vector 
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We thermally perturb the system by a harmonically vary- 
ing heat current density jge'^' with s ~ ico at the surface of 
radius ri . If we impose the boundary conditions of van- 
ishing fields at infinity and a hard core, u{ri) =0 then 
the coupled solution is 



u{r) = 
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The total thermal impedance thus becomes 
^_ 5T{ri) 1 1 
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It should be noted that in solving the same problem on 
the basis of the ordinary heat diffusion equation (|9]l one 
arrives at ( fT9] l but with Cp instead of c;. It is thus seen 
that in doing specific heat spectroscopy by effusion in 
a spherical geometry one obtains again the longitudinal 
specific heat and not the isobaric specific heat. 

One can also consider the case of a soft core, a,., {r]_ ) = 
0. Although the displacement field is altered compared 
to the case of a hard core, the expression for the thermal 
impedance is still found to be given by ( fT9] l. 

The DC-limit gives the heat conductivity. 
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The high-frequency limit is in concordance with the one- 
dimensional result. 
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since short thermal waves cannot "see" the curvature of 
the sphere. It is seen that effusion in spherical geome- 
try in fact gives information on two properties, the heat 
conductivity and the heat capacity, whereas the unidi- 
rectional effusion only gives the effusivity. This is be- 
cause a characteristic length scale, the radius of the heat- 
producing spherical surface, is involved. Effusivity in 
spherical geometry thus makes it possible to derive the 
heat capacity absolutely. However the practical usable 
frequency range will be more limited for a given sen- 
sitivity since the contribution from c/ in fT% will van- 
ish at low frequencies. At high frequency the possibility 
of modelling the contribution to the thermal impedance 
from the heat-producing device itself will also put a limit. 

In real plane-plate effusion experiments the finite 
width of the plate gives rise to boundary effects when the 
heat diffusion length becomes comparable to the plate 
width. The deviation from the simple formula (O is de- 
pendent on the ratio between these to quantities. Since a 
length now appears in the problem this deviation again 
gives the possibility of determining X separately. This 
has been addressed perturbatively |7] on the basis of the 
ordinary heat diffusion equation but not with the 
more exact coupled thermomechanical equations ^ and 
(I?). Thus in fact it seems that of the two simple idealized 
models - the planar and the spherical - of heat effusion 
including the thermomechanical coupling the spherical 
may be the one that mostly resembles its practical reali- 
sation. 



CONCLUSION 

These examples - the unidirectional and the spherical 
geometry - seem to show that it is inherently difficult 
to get the isobaric specific heat directly from effusivity 
measurements. However another well-defined quantity, 
the longitudinal specific heat can be found. 
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